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We compute the free energy for two rows of localized adsorption sites embed- 
ded in a two dimensional one-component plasma with neutralizing background 
density p. The interaction energy between the adsorption sites is repulsive. We 
also compute the average occupation number of the adsorption sites and com- 
pare it to the result for a single row of sites. The exact result indicates that 
the discretization does not induce charge asymmetry and no attractive forces 
occur. 
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The subtleties of electrostatics in condensed matter theory represent a formidable 
and never ending challenge. One topic of much recent activity, has been the attrac- 
tion between two macromolecules of the same charge [Hj. One mechanism that has 
been proposed invokes charge asymmetry related to the formation of lattices or 
Wigner crystals [I21E]- One problem with this picture is that it will create a dipole 
that is inconsistent with the perfect screening sum rules (Blum et al.[ 4J). While 
the formation of Wigner crystals under special conditions is an experimental fact, 
the question of the large asymmetry in the charge distribution needs clarification. 
As has happened in the past the exact solution of the two dimensional Jancovici 
model can provide an unambiguous answer to the puzzle. The interaction be- 
tween two equally charged lines ( which are charged surfaces in 2 dimensions!) has 
been discussed using the exact solution, and is always repulsive [ 6]. Here, we study 
a discretized version of this problem, namely two lines of discrete adsorption sites, 
where the adsorption potential is given by the Baxter [|7] sticky potential. To do this 
we extend the localized adsorption model of a single line to the case of two lines 
of discrete adsorption sites. This extension is non-trivial, and as in a similar case 
discussed in the past, has a simple solution for what we would call a 'commensurate' 
lattice [HI], namely the spacing of the adsorption sites is such that the background 
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charge of the enclosed area corresponds to an entire number of discrete charges. 



1. Formalism 

1.1. Modeling the Adsorption 

Following Rosinberg et al. [ |H| , the adsorption potential for a sticky site located 
at the origin, given by u a (r), is modeled as 

exp[-/K(r)] = l + Ar5(r), (1) 
where A is a positive constant that measures the strength of the adsorption potential [ 

m 

The partition function for a system of adsorption sites with locations given by 
the vectors R m is given by 



' J i=i 



M 



m=l 



d\ , (2) 



where Vq is the potential energy of the one-component plasma in the absence of 
adsorption sites. Expanding in powers of A, it has been shown [|H] that the partition 
function can be written in terms of the n-point correlation functions as 



Z N^2 ^2 ~jP( R rn 1 , ■ ■ "R-rnJ, (3) 



Si 

ml,m2,--m.=l 



where Z% is the partition function of the unperturbed system. 

The difference in free energy from the unperturbed system is the logarithm of 
Z N /Z^, and is given by 

AF = -k B Tj2^jT s , (4) 

«>i s - 

where 

M 

T s = p T (R mi , • • ■ , R m2 ), (5) 

mi,---,m s =l 

and where pt gives the truncated n-body correlation functions, 

pr(ri) = p(ri) 
p T (ri,r 2 ) = p(ri,r 2 ) - p(ri)p(r 2 ) (6) 
p T (r u r 2 ,r 3 ) = p(r 1; r 2 , r 3 ) - p(r 1; r 2 )p(r 3 ) - 
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Figure 1. Two infinite lines of adsorption sites separated a distance d and having 
periodicity a. 



1.2. Correlation Functions of the One-Component Plasma 

The exact solution to the one-component plasma found by Jancovici [ |3] for 
coupling parameter T = Z 2 e 2 /(fcgT) = 2, where Z is the ion valence and e the 
elementary charge of an electron. The n-point density correlation functions are given 
by 

p(r 1; •• -r n ) = p n det 

Here, p is the background charge density, z = x + iy where x and y describe coordi- 
nates on the plane, and z its complex conjugate. After some algebraic manipulation, 
we can rewrite this expression as 

n 

p(r h ■ ■ -r n ) = p n sg 11 ^) 11^°^' " r °U))Px( T j' r <r(i))' ( 8 ) 

O-GSn j=i 

where S n is the group of permutations on n letters, 

Po(r)=exp(-^r 2 ), (9) 

and 

Px( r , r ') = exp [— mp (r x r') • z] . (10) 

The truncated correlation functions, p T (r 1; • ■ • , r n ) are computed by restricting the 
sum over a in Eq. (JEJ) to n-cycles. 

We suppose that the two lines of sticky sites with periodicity a are separated by 
a distance d (see Fig.QJ. We describe the sticky site locations by introducing integer 
variables n» and Ising variables 5i that take on a value of either or 1. Then any R» 
can be written as 

Ri = ariix + dSiy. (11) 

Since we are calculating the sums over all the positions of the particles, all n- 
cycles are equivalent by a suitable relabeling of the summation indices. This leads 
to the general expression 

T s = p s ( s -l)!(-l) s ^^po(R mi , 5l -R m2 , 52 )po(R m2 , 52 -R m3 ,5 3 )---(12) 

ni,— Si,— 

x Po(R-m s ,(5 s — Rmi,5 1 )Px(Rmi,5 1 , Rm2,<5 2 )Px(Rm2,(52 5 Rm 3 ,i 3 ) ' ' ' 
Rm s ,8s ; Rmi ,5i ) • 



-irp(\z tl \ 2 + \z 1 \ 2 )/2+Tvpz,_ t z 1 



lM,7=l, 



(7) 
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2. Free energy 

We first consider px(Ri, R2), and substitute in the adsorption site positions from 
Eq. (JTTJ). This gives 

Px(Ri, R 2 ) = e -™ pad{niS2 - 5im) . (13) 
For the particular choice of background charge density 

p = 2m/ (ad) (14) 

where m is a positive integer, px = 1 when evaluated on the adsorption sites. We 
therefore specialize to densities where this simplification occurs. 

We also find that (Ri — R 2 ) 2 = a 2 (n x — n 2 ) 2 + d 2 (5 1 — 5 2 ) 2 , leading to 

T s = p s ( s - l)!(-l) s ^e^ [(ni -" 2)2+ - +(ns " ni)2] 

ni,— 

e -*'^i-52) 2 + -+(5 s -5i) 2 ] ; ( 15 ) 

where t = irpa 2 /2 and t' = 7rpd 2 /2. Since the sum over the and the sum over the 
Si decompose, we can evaluate the sum over rj, using a transfer matrix. We define 
the transfer matrix, T, to have the components 

T = ^ e -^ 2 /2 ! J • ( 16 ) 

Then 

T s = p s (s - l)!(-l) s tr(T s ) e-* [(ni - n2)2+ - +(ns - ni)2] . (17) 

m,— 

Diagonalizing T gives the eigenvalues X± = l±e~ npd2 / 2 , allowing us to take the trace 
easily. Notice that the decoupling of the Ising variables, rjj and integer variables, rij 
only decouple at densities given by Eq. ()14|). At other densities, these additional 
couplings between the rii and 5i complicates the evaluation of the transfer matrix 
trace. 

The sum over the rii can be expressed in terms of Jacobi theta functions [ 8J, 
where the Jacobi theta function is defined as 

00 

fls(C,t)= E e-*" 2 e 2 ^. (18) 

71= — CO 

First, notice that 

j d( 9 s 3 ((,t)= e -* [n2+ ^ + - +(ni+rt2+ - )2] (19) 

for s > 1 and f d( 3 ((,t) = 1. This leads to the expression 

T s = p s (s - l)!(-l) s [(1 + e-^ 2/2 ) s + (1 - e-^ 2/2 ) s ] / d( 9 3 (C, 7rpa 2 /2) s , (20) 

L J Jo 
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Figure 2. A typical set of free energies Af — A/d_>oo with background density 
pado/2 = 1 (circles, solid lines) and pado/2 = 5 (stars, dotted lines). Here, d is 
measured in units of a and is an integer multiple of 0.1a. We further set Aa 2 = 100. 
The lines are guides to the eye. 




Figure 3. A typical set of free energies Af — A/^oo with background density 
pado/2 = 1 for values of Aa 2 = 0.001 (squares, dashed line), 0.1 (stars, dotted 
line), 10 (circles, solid line). Here, d is measured in units of a and is an integer 
multiple of 0.1a. The lines are guides to the eye. 



Substituting this into equation (jlj) leads to a sum of two series, both of which are 
absolutely summable when \Xp(l + e~ 7Tpd2 ^ 2 )9(0,t)\ < 1. By analytic continuation, we 
extend this sum to the full range of parameters, leading to the free energy difference 
between the OCP with and without adsorption sites given by 



A/ 



knT 



d( In 



1 + Ape- ( ^ 2/2) # 3 (C, 7rpa 2 /2)| - AVe _7rpd2 ^(C, 7rpa 2 /2) 



(21) 



In the limit that d —>■ oo, we expect the free energy to be a sum of the free 
energies of two independent lines of sticky sites. Indeed, this limit yields 

/>/i/ fl d( ln[l + A^ 3 (C,vrpa 2 /2)]. (22) 



A/, 



d— »oo 
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Figure 4. The fraction of occupied sites as a function of Aa 2 for two lines separated 
by a distance d = 0.1a with pad/ 2 = 1 (solid). This is compared to the fraction of 
occupied sites for a single line of adsorption sites at the same background density 
(dashed line). 



In the opposite limit, d — > 0, we expect the free energy to agree with that of a single 
line of adsorption sites with a potential given by 2A. It is easy to see that the free 
energy in this limit is 

A/^o = ~- f 1 d( In [1 + 2\p9 3 ((, npa 2 /2)} . (23) 
a Jo 

Written in terms of the dimensionless constants t = 7rpa 2 /2 = na/d, t' = 
irpd 2 /2 = ird/a and A = Xp = 2X/(ad), we find the change in free energy as the 
adsorption sites approach each other, AF = Af — A/<2->oo, to be given by 

AF = 1» 
« 2 Jo 



l + A(l-e-*> 3 (C,t) ' 
1 + A0 3 (£*) 

Recall that this free energy is only valid when pad/2 = m for any positive integer 
m. Thus, we can compare the free energy of states with the same lattice constant 
a and background charge density p, only for integer multiples of some specific valid 
separation d. To be more specific, suppose we have the free energy at some density 
such that pad/2 = 1 and separation d, then at separation md in the same background 
density, we have pad/ 2 = m. Thus, the free energy formula Eq. ()21|) will be valid 
only for integer multiples of d. 

In Fig. |2l we plot AF as a function of d for two values of the background density, 
p, given by pad^/2 = 1 and pad^/2 = 5, where do = 0.1a is the smallest value of 
d. The lines of adsorption sites are always repulsive. For larger d, AF becomes zero 
quickly. In Fig. El we plot AF as a function of d for different values of Aa 2 . As the 
depth of the potential is increased, the repulsive strength of the interaction is also 



l + A(l + e - t )g 3 (C,t) 
l + A0 3 (C,t) 



(24) 



rp pi 

T- / d( In 

Jo 
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increased. This is indicative of the adsorption sites pinning charge on them, leading 
to their repulsion. 

Finally, we compute the average occupation number of a site. This is given by [ 

IE! 



2, d Af 



{n) = - a X dXk^f- (25) 

This is plotted as a function of Act 2 for in Fig.0]for two lines with separation d = 0.1a 
and pad/ 2 = 1 (solid line) and for a single line of adsorption sites (dashed line). 

It is clear from Fig. 0] that the repulsion inhibits the adsorption of the ions. 
However, the separation d — 0.1a is very small. At separations on the order of the 
site spacing, there is no appreciable difference in the fraction of occupied sites as a 
function of A. 



3. Average density 

In this section, we will compute the average counterion density at an arbitrary 
point Rn. We can find the density directly by fixing the position of one of the ions. 
This will require the computation of the quantity 

s 

T S (R ) = p s+1 ]T II P°( R ^ -**V))Px{B.j, R CT(i )), (26) 

Ri R s J=° 

where the sum is over all s-cycles a. Then the average density is given by [|H] 

°° A s 

(p(x ,y ))-p = J2- T s(' R o) (27) 

where x and y are the components of R . 

For irpad = 27rn, n an integer, and for Rj = anjx + d5jj), equation ()26|) decom- 
poses into the product 

T s = p s+1 s\(-l) s T^T^, (28) 



where 



rp(X) _ e ~7T P x^/2 e ~{npa 2 /2)[nl + (n 1 ~n 2 ) 2 +--+(n s - 1 -n s ) 2 +n 2 s ] 

m—n 3 

y^ e -npaxo(ni+n 3 )+iTTpayo(ni-n 3 ) (29) 



and 



T s (2) = V e -^pyl/-2 e -^ P d 2 /2)[s 2 +{&i-&2?+-+s 2 s ] 

x e~ npdyo ( Sl +Ss ^ e in P d ( Sl ~ Ss ) (30) 
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Figure 5. (p) — p as a function of i/q for xq = in units of a. Here pad = 2n and 
A/9 = 8. We show d = 0.25a (solid), d = 0.5a (dashed), and d = la (dotted). 



Defining z = xq + iyo and 2 its complex conjugate, we find that 

rp{l) _ e -TTpxl/2 e -t{nl+---+n2+(n 1 + ---n s ) 2 ] e -TTpa(n 1 z-n s z) (g-g 

ni---n s 

Using the transfer matrix T and defining a new diagonal matrix, 

= ( n „-TTpd 2 /2+irpdyn+iiTpdxn J (^2) 



Q g-irpd 2 /2+npdyo+iwpdxo J ' 

and its complex conjugate .M, we find 

T j2) = e --p^/2 tr ( rs -i_^ T ^) (33) 

for s > 1 and T (2) = e - npy ^ /2 . Evaluating T s (2) gives 

= [x*-iG + (z) + yj'G^z)] , (34) 

where 4G + (^) = exp(-7rp?/ 2 /2)(A + |l+A| 2 +A_|l-A| 2 ) and4G_(z) = exp(-7rp?/ 2 /2)(A_|U 
A| 2 + A_|_ 1 1 — A\ 2 ), and A = exp(— npd 2 /2 + 7fpdy + mpdxo). 

Using the results of Rosinberg et al. [ IB], we can compute T s in terms of 
F(C, z, t) = {6 3 [C+izt/(na), t]6 3 [(+izt/(ira), t]+6 3 [(-izt/(na), t]6 3 [C-izt/(na), t}}/2. 
This gives 

T s = (-l) s s\p s+1 [X^G^z) + X^G^z)] eT^ 2 fd( 9 s 3 -\(,t)F((,z,t). (35) 

Jo 

Using equation (J2ZJ), we find the average density 

(p(z))-p = -XG + (z)p 2 e-^ f\ F(x -'- f) 



l + \p\+6 3 ((,t) 



., G _ (Z)P ^J\ T J^ J) . (36, 
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Finally, we note that 

AG + {x , y ) = A + j e -^o + e --p(s/o-rf) 2 + 2e -My?Myo-d) 2 ]/2 cos ( vrpc / ;ro )| (37) 
+A_ je-^o 2 + e -p(yo-d) 2 _ 2e -n P[ yi+[y Q -df]/2 cos ( npdxQ ^ ; 

and 

4C7_(xn, y ) = A_ je" 7 ^ + e^^ ^ 2 + 2e - 7r ^o+(w-«0 2 ]/2 cos ( vrp rf Xo )| (38) 
+A+ je" 7 ^ + e ^p(y -d) 2 _ 2e -My 2 +(y -d) 2 ]/2 cos( ^ o) } _ 

Notice that, for npad = 2im, the periodicity of G±(xo,yo) in Xq is always com- 
mensurate with the lattice spacing a. In Fig. |31 we plot (p) — p as a function of yo 
for xo = 0, pad = 2n, and Ap = 8 and for a variety of different spacings. Notice that 
the density has two peaks for large separations but a single peak as the separation 
d becomes smaller than the periodicity a. As a function of Xq, the density is always 
periodic with period a. 

4. Discussion 

The main conclusion of our calculation is that for the geometry that we have 
chosen no attractive forces are induced by the discrete structure of the charged line. 

The derivation of AF is valid for d taking values that are integer multiples of 
2/(ap). Further, Eq. ()21|) gives the correct free energy for d = (as we have already 
seen). It is conceivable, then, that Eq. (|2~Tj) is correct for all values of p and d. 

One of the features of this model that makes an exact evaluation possible is that 
a density p can be found such that the Ising variables Si and the integer variables 
Hi are uncoupled. When pad 7^ 2m for an integer m, a coupling does indeed arise 
that makes the computation of the free energy more difficult. Additionally, if the 
adsorption sites are not aligned, an additional x Si Aa component arises in Rj, where 
Aa measures the degree of misalignment. This component introduces a coupling 
between the rij and Si that will be discussed in future work. 
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